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Conical brane singularities in six-dimensional flux compactification models can be resolved by 
introducing cylindrical codimension-one branes with regular caps instead of 3-branes (a la Kaluza- 
Klein braneworlds with fluxes). In this paper, we consider such a regularized braneworld with 
axial symmetry in six-dimensional Einstein-Maxwell theory. We derive a low energy effective theory 
on the regularized brane by employing the gradient expansion approach, and show that standard 
four- dimensional Einstein gravity is recovered at low energies. Our effective equations extend to the 
nonlinear gravity regime, implying that conventional cosmology can be reproduced in the regularized 
' braneworld. 
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I. INTRODUCTION 



Recent developments in string theory provide us with a remarkably new picture of our universe, in which our universe 
is described by the motion of a 3-brane in a higher dimensional spacetime. In addition, it has been suggested that 
I • extra dimensions warped due to self-gravity of branes may play a key role in addressing the gauge hierarchy problem. 
Q-i' A simple toy model of warped braneworlds was proposed by Randall and Sundrum [l[ and its various consequences 
^ • has been studied extensively Qj- The idea of warped extra dimensions is plugged in warped flux compactiflcations, in 
which reliable inflation models based on string theory has been constructed [3| (see Ref. Q for a review). 

While the flvc-dimcnsional (5D) Randall-Sundrum-type models [J are simple enough, offering us a feasible setup 
for studying the behavior of gravity in braneworlds, brane models in six or higher dimensions have a much richer 
^ , structure such as flux-stabilized compactiflcations. In particular, codimension-two braneworlds with football-shaped 
■ extra dimensions attracted much attention initially, due to their potential mechanism for solving the cosniological 
ly-j \ constant problem [1, [1,01- However, a higher codimension brane develops a conical or worse singularity once taking 
' into account its self-gravity, which hampers a substantial analysis of higher dimensional braneworlds. Thus, it is 
necessary to regularize higher codimension branes in a well-controllable way to construct reasonable models of higher 
dimensional braneworlds. One of the regularization procedure is to replace a higher codimension defect with a 
codimension-one brane, whose extra dimensions are conipactified in the traditional Kaluza-Klein manner. This was 
, done in an specific model of 6D flux compactiflcations in d, [^, [l^ . This type of models may be regarded as a variation 
L; ■ of Kaluza-Klein /hybrid braneworlds [H IH, [H, [11, [H . Only pure tensiom gravitational Shockwaves, and black holes 
. ^ , have so far been known to be accommodated on codimension-two branes [161] , but with this regularization scheme one 
can put arbitrary energy-momentum tensor on the branes. Performing perturbation analyses, it has been shown that 
5h weak gravity sourced by arbitrary matter on the regularized brane reproduces standard 4D one 

In this paper we further investigate the behavior of gravity on the regularized brane in the 6D flux compactiflcation 
models, generalizing the previous result to cover nonlinear gravity. We keep in mind the bulk-brane conflguration 
of [1, but start with a more general metric ansatz. Employing the gradient expansion approach, we solve the 
6D Einstein-Maxwell equations to derive an effective theory on the brane. The gradient expansion approach in 



braneworlds is based on the (4 -I- l)-decomposition of the bulk spacetime [18[ and was originally developed in the 
context of the 5D Randall- Sundrum model p^. [20j. The gradient expansion has been widely used in various brane 
models (e.g., [2l[). Recently, along the direction of the regularized braneworlds, the geometrical projection approach 
of [isj has been generalized to non-Z2-symmetric models in an arbitrary number of dimensions 

The plan of the paper is as follows. In the next section we perform the (5 -f l)-decomposition of 6D spacetime. 
In Sec. mil the evolution equations along one of the extra direction are solved by invoking the gradient expansion 
approach up to flrst order. We explain how conical singularities in the bulk are regularized by introducing extended 
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4-branes in Sec. IIVI In Sec. |V] we derive the effective tlieory on the regularized branes using the result of Sec. IIIIl In 
Sec. IVII wc draw our conclusions. 

II. BASIC EQUATIONS 

The 6D action we consider is 



S ^ / 



(1) 



where = FmnF^^^ and Fmn = dnAN — OnAm is the field strength of the U{1) gauge field Am- M is the 
fundamental scale of gravity and l/2Lj is the cosmological constant. Branes will be added in Sec. |TVf and for the 
time being we will concentrate on solving the 6D geometry. We assume the axisymmetric metric ansatz in the form 
of 

gMNdx^'dx^ = L]^ + £2e2V'(-.s/)/(y)rf02 ^ (x, y)dedx^' + y)dxf^dx-', (2) 

f{y) 

with 

9tiv = o?{y)ht,^{x,y). (3) 

The capital Latin indices numerate the 6D coordinates, while the Greek indices are restricted to the 4D coordinates. 

We write the 6D Einstein equations in the form of the evolution equations along the y-direction and the constraint 
equations. The evolution equations are 

nyO.K.^ + KRy ^ 'R^ - -^Sy - ^ {f^mF'^' ~ \5yF^ , (4) 

where = ^/J /Li, Ky is the extrinsic curvature of ?/ = constant hypersurfaces, and K is its 5D trace. ^Ry is the 
5D Ricci tensor. Here, fi ~ n and 9. The Hamiltonian constraint is 

'R + K^Kf - = -1 - [FyMPy^^^ - \f^^ . (5) 

and the momentum constraints are 

'D, [xy - syx) = -l^F,,MFy'^'ny, (6) 

where ^Dj^^ is the covariant derivative with respect to the 5D metric and Uy ^ Lj/ ^f] . 
The Maxwell equations are given by 

Vj^F^^^ = 0, (7) 
where Vat is the covariant derivative with respect to the 6D metric. 

III. GRADIENT EXPANSION APPROACH 

We employ the gradient expansion approach to solve the 6D bulk geometry [l^ [l^]- Assuming that I is not so 
different from L/, the small expansion parameter here is the ratio of the bulk curvature scale to the 4D intrinsic 
curvature scale, 

e^i^\R\. 



^ There we will glue the solutions with different cosmological constants together, and so the subscript I will eventually indicate the 
different regions of the spacetime. 
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The metric and extrinsic curvature are expanded as 

h - ;,(o) + 4- . . . 



(0) (1) 



(0) (1) 



We may assume naturally that 



6. = £^/'6fy') + --- 



1 /o <^^^' 



(8) 



The (y^) component of the field strength is expanded as 



(0) (1) 



(1/2) 



Since dfj^Ag ~ e^^^dyAg, we have -F^e = E^^^'^F^e + ■ • ■ • We also assume that 



(1/2) 



(1) 



(9) 



Therefore, the F^aF'^" part of the bulk energy- momentum tensor will be higher order in our treatment and it does 
not contribute to the low energy effective theory. The 5D Ricci tensor is given by 



5t> e 



(10) 
(11) 



and ^R^ — ©(e^/^),^ where stands for the 4D Ricci tensor of the zeroth order metric and I?^ is the 

covariant derivative with respect to . 

A. Zeroth order equations 

At zeroth order in the gradient expansion, the evolution equations and the Hamiltonian constraint are 



Li 



(0) / 



(0) (0) (0) 

K,' + Kg') k; 



(0) (0) (0) 



(0) (0) / <"> \2 (0) (0) 



2 K^'K^ 



1 



3 



ye 



S " 



(e-^'"' F,^ 

(0) s 2 



(0) s 2 



— (e-^'"' Fg] 



(12) 
(13) 
(14) 



The Maxwell equations at zeroth order is 



%{^^a'e-^'"'FZ)^0. 



(15) 



The above equations are solved by 



a{y) = y, 
/(y) = -:^ + 



20 2/3 i2y6 ' 



(16) 
(17) 



^ It turns out in Appendix Rl that we in fact have ^Rf' = ©(e^/^). 
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and 

Fyg=£M'^^e^'"\ h^^}(x,y)^h^,{x), ^t"' 2/) = V'^"^^), (18) 

with and Q being constants that characterize the geometry of the internal space. 

We assume that the metric function f(y) has two positive roots, y^ and ys{< yN)- In other words, we choose the 
parameters fi and Q so that f{y) has two positive roots. We impose that Aq vanishes at these poles, yielding (in the 
northern half) 



Vn 



From this we obtain F^g= df^ip^^^ • A^g \ Substituting this into the /i component of the momentum constraints, we 
have d^ip''^'^ = 0. Therefore, must be constant and without loss of generality we may set = 0. Furthermore, 



(1/2) 



this implies F^g— and hence F^g — 0(£^/^) + • • • . 

The above bulk geometry is our "background," which is essentially the double Wick rotated Reissner-Nordstrom 
solution [Til,!!!. 



B. First order equations 



Going to first order in the gradient expansion we will obtain the gravitational field equations that govern the 

/QN (l/2)(l/2) (l/2)(l/2) 

behavior of the 4D metric hi^l . The first order equations may contain terms like ^F^yF'^*' and r^KgR^. However, 
in Appendix [XI we analyze the ©(e^/^) equations and show that the gradient expansion of F^j, and Kg in fact begins 
with C'(e'^/^). Therefore, in the following we will drop the contributions from such terms. 
At first order, the {jiv) component of the evolution equations is given by 



Li 



(1) 


(^4 M\ 




\v V ) 



(1) 

k:, 



(1) 
if/ 



Iff- 



(20) 



where we defined the useful combination 



1 / (0) (1) (1) (0) 
^ Jf.yFygFy'' + FygFy' 



(21) 



The 4D Ricci tensor i?^ does not depend on y because it is computed from which is a function of only a;''. The 
traceless part of Eq. ([^0]) is found to be 



(22) 



(1) 



where we defined the traceless part of the relevant tensors as K^'' ■■=K^ -(l/A)5^ and R^'' i?^'' - (1/4)5^'' i?. 
Eq. ([22]) can be integrated to give 



3yV7 ^ 



c;(a:). 



(23) 



where the traceless tensor C^(x) is the integration "constant" to be fixed by the boundary conditions. 
The trace part of the evolution equations is 



Li 



2/ 



A (1) 

y 



(24) 



and the 



component of the evolution equations is 



77 

Li 



dyKg" 



f 



(1) 



2/ 



7^. 



(25) 
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The Hamiltonian constraint at first order reduces to 



y 



1 \/7 



3_ 



(1)^ 4 (1) 



K 



(26) 



It is possible to solve the set of equations p4 | - (|26|) and find the bulk profile of the extrinsic curvature and T (see 
Appendix [B]) . Then, with the help of the relation 



(1) , (0) (1) , 



one finds the bulk profile of the gauge field. Equivalently, one can use the first order Maxwell equations, 



(27) 



(28) 



to solve for the gauge field. As will be seen, however, we do not need the explicit form of the solution for the extrinsic 
curvature and gauge field in order to derive an effective theory on the brane; it suffices to solve the bulk evolution of 
the combination of first order variables 



o (1) (1) 



Li \2f yj 



The evolution equation for K, can be obtained as follows. Using the relation ([27]), it is straightforward to show 



Then, using the evolution equations 



and the Hamiltonian constraint (|26p . we arrive at 



2 

dy (y^^/fJc) = ^LiR. 



The structure of Eq. ([?T|) is identical to the traceless equations , and the general solution is given by 



12^/^/7 



LiR 



(29) 



(30) 



(31) 



(32) 



where x(x) is an integration "constant" to be determined by the boundary conditions. 
The fj, component of the momentum constraints at first order is given by 



^ dyf 



Li \y 2/ 



(1) 1 (1) C' Lr 

V.i)^ —rV^Ae Fy^ -4, 



which can be simplified using K, to 



(33) 



Eq. ([551) the Bianchi identity, V^, \R^ — {\I2)5^ Pi\ — 0, imply the constraint for the integration constants: 



(34) 



IV. REPLACING CONICAL BRANES WITH REGULARIZED 4-BRANES 



The metric function /(y) vanishes at j/at and ^5. These points develop conical singularities in general and they 
are regarded as source 3-branes. To resolve the singularities, we replace each of the conical branes with a cylindrical 
codimension-one brane and fill in the interior with a regular cap . The geometries of the two caps and central 
bulk region are described by the 6D solutions obtained in the previous section with different cosmological constants 
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(L+ (i-) for the north (south) cap and Lq for the central bulk). Near the pole y = yp {p = N, S), we have 
f{y) ~ dyf\y=y^{y — yp). In order for the cap to close regularly at yp, we impose 



e \dyf\ . 

Li 2 



2n. 



(35) 



Clearly, it is required that ^{x,yp) = constant. Without loss of generality we can set this constant contribution to be 
zero. In particular, we have ^^^^ ^ {y ~ Up) near y = yp- 
We consider the following action for each 4-brane 



(36) 



where A± is the brane tension, v± is the vacuum expectation value of the brane Higgs field and Sj- is its Goldstone 
mode. S*^; represents the matter action on the brane. The brane action necessarily couples to Af^ in order to account 
for the jump of the Maxwell field at the brane. This implies that the brane location is fixed[see Eq. p9|) below]. The 
total brane energy-momentum tensor derived from the above action is 



T 



±0 

fi(tot) 



(a^s± - e^^) {d^^± 



eA')~-{d;j:±-cA;) 



d^Y.+ - eA' 



±0 



(37) 



0. 



where T^" is the energy-momentum tensor of the matter fields on the branes. We assume that T^^ = T^'^ 

Each of the caps and the central bulk spacetime are glued together so as to satisfy the Israel conditions [l^l and 
the jump conditions for the Maxwell field. From now on we will suppress the index ± if it is not necessary. We start 
with assuming that the brane location is given by a x-dependent function y = (p{x). The brane induced metric is 

given by qfipdx^dx^ = P f {ip{x))d9'^ + ifP' {x)h^l^}, dx^^ dx^ + 0{e). The equation of motion for each scalar field E± at 
zeroth order reduces to dgTi^^ = 0, and hence 



(38) 



Here, n± must be integer because I]± is the phase of the Higgs field and so e»5^(^'+2T,2;) _ ^iT.(e,x) ^ This property 
should hold at any order in the gradient expansion. Therefore, the solution for E including higher order corrections 

should be of the form ^ nO + a^°'^ {x) -I- ecr^^) (a;) H . 

We shall show now that the brane location is in fact independent of x. To do so, let us consider the 6 component 
of the Maxwell jump conditions at zeroth order: 



(0) ■ 



-ev 



eA 



(0) 



(39) 

iny,0). The 

right hand side is independent of x and hence the brane location must be y = constant {—: y±). The Israel conditions 
are given by 



where [[^4]] :— A\y^ip+^ — A\y=^-^ and n is the unit normal to the brane. At this order, we have 



-M 



yif^^dyi 



Li \y 



2/ 
Li y 



Af4 2A/4 (^^^ 



(0) 



eA 



(0) 



A 



eA 



(0) 



T 



fi{tot) 

e{o) 



B(tot)- 



(40) 
(41) 



The above conditions ([55|) - (|¥T|) determine the brane location y± and the parameter of the solution n±. For the current 
purpose we need no further detail; see Sec. II of Ref. 17] for more explanation of the configuration of the branes. 

V. EFFECTIVE THEORY ON A REGULARIZED BRANE 

We go on to specifying the first order boundary conditions at the poles and branes. As to the regularity conditions 
at the poles, it is required that 



(42) 
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With this, the evolution equations for the extrinsic curvature (|22l) . p4l) . and ([25]) are regular at the poles. We also 

require that | Fyg | < oo at y = y^. This implies that A''^^ ^ (jJ ^ Up) near the pole. Noting that ■f/'^^) ^ iv ~ Up), we 
have /C < |y — yp\^^^ near y = yp. 

The (iiiy) component of the Israel conditions at first order are given by 



(1) 



T,: 



1 



Af4 M4 £2/ 

In deriving this we used {d'^T, — eA'^Y^^^'' — which is shown in Appendix [XI The traceless part of Eq. reads 



(43) 



where is the traceless part of the energy-momentum tensor. Noting that'^ 



(44) 



1 rr ■ 

m]^-[[3K,^+4Kj> 



1 1 



(0) 



ye 



(1) 



V7 W _i 



and using the zeroth order Israel conditions and jump conditions for the Maxwell field, we find that the trace part of 
Eq. (031) yields 



[[IC] 



1 ji A 
4M4 A • 



(45) 



The junction conditions (|44p and (|45p together with the momentum constraint ([33]) imply the local conservation law 
for the energy-momentum tensor of brane localized matter: 



V^T" = 0. 



(46) 



Now we are to fix the integration constants. From the regularity condition at the north pole one can determine the 
integration constant in the north cap and we have 



y^ - Vn 



^ 3y4V7 

Then, using the Israel conditions at the brane we obtain 



i+R; (y+<y<yjv). 



1 



1 



where 



(.yff-y+)L+ + y+Lo^^ 



r iy~ <y < y+), 



y+Vf+- 



and /+ :— f{y+). The regularity at the south pole determines the extrinsic curvature in the south cap as 



{ys<y<y-). 



The Israel conditions at y = y_ imply 

{y'L-yDL^-ylLo 



y-Vf- 



(47) 



(48) 



(49) 



(50) 



(51) 



The gauge field A^^^ is continuous across the brane in order for the brane action to be well-defined. The continuity of the induced 
metric imposes that i/"'^^ is continuous across the brane. 
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which is rearranged to give 

L ^T+'^ + yt^T-/) , (52) 



where 



Vn 

Liy^dy (53) 



and /_ :— f{y-)- Obviously, for the trace part we have the same relation as ([52]) with the substitution R^^ R/A and 
1^'^ — (and thus the constraint for the integration constants ((34|) is trivially satisfied). Taking into account 
that the induced metric on the north brane is given by y\h^^, we define the branc Ricci tensor as 7?,^ :— R^ly\. 
Combining the traceless and trace equations, we finally arrive at 

where we defined the 4D gravitational coupling at each brane as 

-2 - ^± (55) 



and the energy-momentum tensor integrated along the 0-direction as := 27r£y/Y±Tjt'^ ■ Eq. shows that 
gravity at low energies is described by general relativity (when matter on the south brane can be neglected.) This 
generalizes the perturbative analysis of [1, to the nonlinear regime. The effect of the extra scalar mode that 
appears in [l^ is higher order in the gradient expansion and hence is not observed in the present leading order 
analysis. 



VI. CONCLUSIONS 



In this paper we have considered six-dimensional Einstein-Maxwell model of warped braneworlds, where the extra 
dimensions are stabilized by a flux. Conical singularities in such models can be regularized by replacing the 3-branes 
with codimension-one branes. Using the gradient expansion approach, we have derived a low energy effective theory 
on the regularized brane, and have shown that standard 4D Einstein gravity is reproduced on the brane. Though the 
present analysis is restricted up to first order, one can go on to the higher order expansions to include the effect of 
Kaluza-Klein modes. The scalar zero mode corresponding to the fluctuation of the internal space volume would also 
appear at second order. Computing higher order corrections would be interesting and is one of the remaining issues. 

Our result has been derived without relying on linear perturbations in the Minkowski brane background. The 
effective Einstein equations (|54p are valid even in the nonlinear regime, and hence can be applied to the cosmological 
dynamics at low energies (i.e., to late-time cosmology in which the Hubble horizon is much larger than the com- 
pactification scale). Since the effective theory coincides with 4D general relativity, the present model can reproduce 
conventional cosmology on the brane. Recently, cosmology on a moving brane in a given (static) background has 
been investigated in the context of regularized braneworlds with flux compactifications [1^, [2^ . Their results are not 
compatible with the realistic cosmological evolution, suggesting that the static bulk assumption is an oversimplifica- 
tion. On the other hand, we have started with a more general metric ansatz ([2]) whose internal space component is 
nontrivially dependent on the external spacetime coordinates. Thus, constructing explicit cosmological solutions in 
the 6D flux compactification model is not so simple as in the 5D Randall-Sundrum-type braneworlds (see, however, 

diii). 

6D braneworlds are often discussed not only in a simplified Einstein-Maxwell setup but also in the context of gauged 
chiral supergravity (29l.[30j. Supergra vity models generally have conical brane singularities too, and the regularization 
procedure has been proposed in |3l|. The present gradient expansion approach is useful for the supergravity 
braneworlds as well. 
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APPENDIX A: Oie'^^^) QUANTITIES 

Starting with the assumptions ([8]) and ([9]), we shall show in this appendix that the leading order terms in the 
gradient expansion of Kg'^ and are in fact 0{e^^'^ ). 
The fj, component of the 0{e^^^) Maxwell equations reads 



/ , (1/2) 



0, 



(Al) 



and thus we have 



(A2) 



The ©(e^/^) evolution equation reduces to 



1 1 



(1/2) 



1 (0) (1/2) 

— — Fo F^y 



(A3) 



which can be integrated to give 



(1/2) 



(A4) 



(l/2)(l/2) 

The 0{e) evolution equations contain terms like F^j,i^''^(x hfj^xC^Ci /[y f{y)]. Thus, in the cap regions it is required 
that Ci = because otherwise this term would show a singular behavior at the poles. Similarly, to ensure the regular 

(1/2) 

behavior of Kg'' at the poles we impose C2 = in the cap regions. To fix the integration constants in the central 
bulk region, we invoke the junction conditions at the branes. The jump conditions of the Maxwell field and Israel 
conditions imply, respectively. 



ny 



(1/2) 



M4 



(1/2) 



^(9.1](°)-e4")) (a-S-eA-) 



(1/2) 



(A5) 
(A6) 



(1/2) (1/2) 



Combining these two equations and noting that Kf'= —F^f^ in each cap, we obtain two linear algebraic equations 
for Cj; and C^: 



■(1/2) 



(1/2) 



0. 



(A7) 



Thus, it is now clear that both and C2 vanish in the bulk. 
To sum up, we have shown in this appendix that 



, (3/2) 

X/ = £3/2 



£3/2 pyfi _|_ ^ 



(A8) 
(A9) 



and that 



(a^E - eA^)^^/^^ = on the branes. 



(AlO) 
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APPENDIX B: SOLVING THE FIRST ORDER EVOLUTION EQUATIONS AND HAMILTONIAN 

CONSTRAINT 



For completeness we present here how we can solve the evolution equations ([24]) . ([25l) . and the Hamiltonian con- 

(1) (1) 

straint (pS)) analytically. Eliminating JF, we obtain the coupled evolution equations for and K^: 

(1) /O fl f\ (1) A (1) 

+ - 0> (Bl) 

(1) 

Then, we substitute Eq. (jBip into Eq. (IB2p to obtain the single second order differential equation for K-;^. In terms 
of 

(1) 
K ^ 

^ — — / ^ , s , (B3) 

the equation can be written in a simple form: 

d^n - §^dyn = 3 ^ LiR, (B4) 

where 



The general solution to Eq. (|B4[) is given by 

n^C,ix)+C2ix)Iiy)+3LiR dy dyliy) % (B6) 

i ^ y^fdy{Vf/y)dyI{y) 

Eqs. (jBip and (|B3[) allow one to write the extrinsic curvature in terms of il. Then J-' can also be written in terms of 
using the Hamiltonian constraint. 
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